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Abstract 

A simple assumption of an emergence in gas of small atomic clusters consisting of c particles 
each, leads to a phase separation (first order transition). It reveals itself by an emergence of 
"forbidden" density range starting at a certain temperature. Defining this latter value as the critical 
temperature predicts existence of an interval with anomalous heat capacity behaviour c p oc AT -1 / C . 
The value c = 13 suggested in literature yields the heat capacity exponent a = 0.077. 
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I. INTRODUCTION 



The theory of gas-liquid condensation is, probably, the most famous unsolved problem 
in the classical statistical mechanics [T]. Numerous attempts to attack the problem have 
been made during the last hundred years. They were based on a wide range of different 
techniques: from cummulant expansion to field theoretical methods of phase transitions [2] . 
A considerable step in this direction had been made by the cluster (droplet) theory of 
M.E. Fisher [3j. This theory predicts an essential singularity of the free energy at the 
condensation point. 

A simple model of condensation which opens the way for appearance of a critical point 
and the corresponding phase separation is suggested here. This model reveals the basic 
desirable features of the condensation and allows a new and self-consistent definition of the 
critical point. Moreover, it identifies the famous heat capacity singularity and explains it up 
to the calculation of the divergency exponent in an excellent accordance with the measured 
data. 

Isolated clusters of atoms and molecules have been observed experimentally in molecular 
beams and studied theoretically pi]. Stability of such clusters has been studied also in a 
liquid-like environment by S. Mossa and G. Tarjus [5]. They have shown that the locally 
preferred structure of the Lennard- Jones liquid is an icosahedron (13 atoms), and that the 
liquid-like environment only slightly reduces the relative stability of it. 

Scattering experiments can also be regarded as an additional indirect argument in favor 
of clustering in liquids. For example, argon radial distribution function [6] shows neither 
temperature nor density dependence of its first maxima abscissae, i.e. internuclear distances 
in solid, liquid and gaseous argon are inherent characteristics of the material. In other words, 
this phase independence can be attributed to the persistence of small dense clusters. 

More detailed study of experimental evidences in favor of the existence of relatively stable 
small atomic clusters will be published elsewhere [7]. 

II. BASIC ASSUMPTION 

Therefore it is possible to formulate the following basic assumption: elementary particles 
of gas (atoms or molecules) form small, relatively stable clusters consisting of c particles 
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should be considered as (at least) two-component system (see Fig. [T|. 

The ground state of the system under consideration is expected to be a full separation as 
the energetically preferable one (we do not address here those special cases when geometry 
allows packing denser that FCC or HCP). On the other hand, at high temperature the 
system remains a mixture of atoms and clusters. Thus, at finite temperature the separation 
into two phases occurs. 

This observation helps us to answer a very natural question: why do we suppose only 
one size clusters to be formed or, at least, to be stable. Unfortunately, we do not know an 
a priori reason for this. On the other hand, as one sees, the existence of clusters of one 
size leads to the separation. Therefore, the existence of clusters of any different number 
of particles would reveal itself through multiple separations. To the best of the author 
knowledge, it is not what happens in Nature with simple liquids. Thus, this a posteriori 
argument justifies our basic assumption. By the way, one may attribute complicated phase 
diagrams of complex liquids to the existence of clusters of different sizes and nature. 

Such a model reveals a universal behaviour. Indeed, a close vicinity of the critical point 
(if it exists) has to be governed by the universal properties of the two-component mixture 
separation, regardless of the specific details of the inter-particle interaction. The latter 
affects the critical parameters, i.e. physical coordinates, but not the system's behaviour. 

Our basic assumption plays a role analogous to the Cooper pairing in the superconduc- 
tivity theories: it is a microscopic phenomenon underlying the macroscopic one. Knowledge 
of the exact (probably, quantum) mechanism of this clustering is not crucial to understand 
the liquid-gas transition. 

III. FREE ENERGY 

We start with the expression for the Helmholtz free energy for a two-component slightly 
non-ideal gas mixture [H] 



Let Ni = n be the number of clusters containing c particles each, N2 = N — cn and N is 
the total number of particles. (3 = (/c^T) -1 as usual. As it said, we assume that all the 
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clusters have the same and constant number of constituent particles, c. Aj = (27i(3/m i ) 1 ^ 2 h 
is a thermal wave length and mj is a particle mass. Eg stands for a cluster binding energy 
and Bij denote second virial coefficients. Thus, 
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and the internal energy 



3 1 



V ^ 7y 2/3 



[AT-(c-l)n] + £ B n+-^(/3; 



where 



(3) 



(4) 



B{/3; n) = B xx {fi)n 2 + 2B 12 ((3)n(N - cn) + B 22 ((3){N - cnf. 
Within the same approximation (slightly non-ideal mixture) the equation of state reads [8] 

(5) 



P0=I[JV-( c -l)n] + ^B(#n). 



A dynamic equilibrium configuration of the two-component system is defined by the value 
of n corresponding to the minimum of the total free energy. Simple differentiation of Eq. (J5]) 
leads to the main equation determining n: 



In (A? 



cln{\* N - Cn 



1 



. 1 ——. + (3E B + v B> n (P;n) = 0, 



(6) 



or 



ln(A 3 xp) - cln (A 3 p(i - cx)) - § lnc + (3E B + pB' x (/3; x) =0, (7) 

where p = N/V, x = n/N, A = A 2 , Ai = c _1 / 2 A and B(/3;x) = B u x 2 + 2B 12 x(l - cx) + 
B 22 (l-cx) 2 . 

One has to solve analytically Eq. Q, i.e. to find x = x(p). Instead, we found an 
inverse function, p = p(x), where x 6 [0, 1/c]. It is easily done with the aid of Lambert 
W- function [H] (another notation: w-function): 



where a = c 3//2 exp^Es). In fact, the equation of state (|5j) in the form 



ax 
'1-cx) 



exp < -W 



ax 



_(l-cx) c \ (c-l)A 3 



(8) 



Pf3 = p[l-(c-l)x]+p 2 B(f3;x) 



(9) 
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and Eq. ^ define P(p) using the parameter x. The most interesting feature of Eq. ([8| is 
the existence of "forbidden" values for p. This behaviour is governed by the sign of the 
derivative B' x (f3;x). Namely, if for a given j3 it remains negative for all permissible values 
of x, then p ranges over the entire positive semi-axis. It is clear from the behaviour of 
Lambert function in the negative range [9]. If the expression changes its sign to positive, 
an equilibrium solution jumps from the Wo branch, continued from positive argument, to 
the W_ i one. Moreover, the positive range of the expression will have another "forbidden" 
region as the absolute value of the Lambert function's negative argument cannot exceed 1/e. 



IV. THE CRITICAL POINT 



The standard definition of a critical point is 

dP\ fd 2 P 



dp J 8 V dp 2 



0. (10) 



However, this definition is not applicable if one expects some singularity to be revealed at 
this point. Moreover, as we just saw, there exists some special behaviour characterized by 
the sign of B' x (f3;x). Thus, the very last (critical) point before the p-axis becomes "teared" 
up is defined by B' x (f3 c \ x c ) = 0. In fact, this equation defines critical parameters: (inverse) 
critical temperature, {3 C , and critical concentration, x c , satisfying 

[B 12 (f3 c ) - cB 22 (/3 c )} + x c [B n (/3 c ) - 2cB 12 (f3 c ) + c 2 B 22 (/3 c )} = 0. (11) 

The left-hand side consists of smooth monotonic functions of (3 (second virial coefficients) 
and is linear in x and, thus, attains its extremum at a limiting point. It cannot be x c = 
because our physical system is supposed to be stable for small concentrations. Therefore, 



the only possibility is x c = 1/c and Eq. (11) becomes 

B n (/3 C ) - cB 12 (p c ) = 0, (12) 
whose root, (3 C , is the inverse critical temperature. Naturally, these equations for x c and (3 C 



are strongly approximation dependent ones. Higher viral expansion will complicate Eq. (11) 
leading to different values for the roots x c and (3 C . 

An important observation to make here is the atom-cluster (B\ 2 ) and the cluster-cluster 
(B 22 ) interactions should be substantialy weak in comparison with the inter-atomic (-Bn) 



one, since part of the gas energy is accumulated in the cluster bindings. It results, in turn, in 
"shallow" potential well with a much shorter repulsive part and relatively small inter-cluster 
distance and, then, in a much higher density of a heavy component of the gas. 

This new definition of the critical point, B' x (f3 c ;x c ) = 0, allows one to write down an 
expansion in the vicinity of this point 

B' x {^x)^B^ c] x c )A^ + B^ c] x^Ax, (13) 

where A/3 = (3 C — (3 and Ax = x c — x. Substituting this, x — > (3 — > (3 C and 1 — cx — > cAx 
into Eq. we obtain the main Eq. (|7| in a close vicinity of the critical point 

c In (Ax) -cA = pB xx Ax, (14) 

where cA = pB x/3 A(3 + (3 c Eb — (c — 1) In (A;?p) + (c— |) In c. This equation is solved as before 
with the aid of the Lambert function and its solution reads: 

Ax = e A exp{-W(-lpB': x e A )} (15) 

with e A = (A^p) 1 c c 1- ^ exp ((3 c Eb + p-B^A/3) }. This looks like an ultimate solution 
of the problem, at least, in the vicinity of the critical point but it does not account for the 
basic feature — the discontinuity of p-scale — and it should be used very carefully. 



V. SPECIFIC HEAT 



The internal energy is given by 

U _ 31 - (c- l)x 
N ~ 2 $ 

and the specific heat — by 



+ E B x + pB'p(P;x) (16) 



p 



= k B {§ [1 - (c- l)x] - p(3 2 B'^} + k B f3{l(c - 1) - PE B - p pB» p } x' p . (17) 

Therefore, if one looks for special behaviour of this quantity in the vicinity of the critical 
point then x and x'p have to be examined. We also make use of the fact that cy on the 
critical isohore behaves like c p in the second order phase transition. [TU] 
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We start with substituting Eq. (13) into Eq. (|8j and note that B" x (f3 c ; x c ) = Bn(/3 C ) — 
2c£> 12 (/3 c ) + c 2 B 2 2((3 c ) = -B 12 (/3 C ) — cB 2 2((3 c ). It represents the cluster-atom and the cluster- 
cluster interactions which are supposed to be very small. Thus, one can expect existence of 
an interval where B' x ((3;x) ~ B x/3 (f3 c ; x c )A/3 and 



In {Xlp) = In 



a c /c 
(cAx) c 



i 

c-1 



W 



a c /c 
{cAx) c 



x B>^{f3 c ;x c 



(c-l)A» 



-A/3 . 



(18) 



Further consideration depends on the sign of B x/3 ((3 c ;x c )A(3. In the homogeneous phase 
it is negative and we are on the Wq branch with a small positive argument. Here it is enough 
to take [9] Wo(y) ~ y and subsequently 



A^P 

The relevant root behaves as 

i 

(cAx) c c - 1 
a c /c 



a c /c 


c-1 


{- 


a c /c 


i 

c-1 


_{cAx) c _ 




_(cAx) c _ 





^B"{p c -x c ) 



A/3 



B'LpSfici x c 



-A/3 or Ax oc (A/3) 1 = . 



(c-l)A3 

It means that the derivative Ax/ A/3 and therefore the specific heat will show here the famous 
dependence c p oc (A/3)~ = . In view of the previous suggestion, c = 13, this exponent becomes 
a w 0.077. 

An analogous calculation cannot be done for a nonhomogeneous phase as an equilibrium 
solution does not exist in this region. 



VI. CONCLUSIONS 



A model that explains basic features of the condensation is presented. A simple as- 
sumption of relative stability of only one type of clusters statistically emerging in the gas 
immediately leads to the first order phase transition (phase separation) at some finite tem- 
perature. It is experimentally observed as a condensation process. 

It should be stressed that this model is by no means a simplified version of Fisher's one. 
As much as the monogamy is not a simplified version of the polygamy and the monotheism 
is not a simplified version of the polytheism. 

Mathematically, the condensation reveals itself as a forbidden density (volume) region. 
The density jumps from its gaseous value to the liquid one. No intermediate values are 
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allowed. A correspondent region for the Van der Waals equation is a well-known S-shaped 
instability. It needs special auxiliary construction to be treated as a metastable state. 

This paper presents a new concept of the critical point: it is a point of the density's 
continuity failure. This definition coincides graphically with the old one but it allows to 
construct a convenient expansion in the close vicinity of the point under consideration. It 
demonstrates the famous singularity with the exponent a ~ 0.077 that agrees excellently 
with known data. 

Acknowledgement 

The author is grateful to Alexander Voronel and Moshe Schwartz for the valuable dis- 
cussions. A financial support of Alexander Voronel during a part of this study is kindly 
acknowledged. Extensive editorial efforts of Ely Klepfish made this manuscript readable. 



[1] Isihara A 1971 Statistical Physics (New York: Academic Press) 
[2] Langer J S 2000 Ann. Phys. 281 941-90 
[3] Fisher M E 1967 Physics 3 255-83 

[4] Haberland H 1994 Rare Gas Clusters Clusters of Atoms and Molecules Chem. Phys. 52 
(Springer Verlag) p 374-95; Berry R S 1999 Phases and Phase Changes of Small Systems 
Theory of Atomic and Molecular Clusters (Springer Verlag) p 1-25 

[5] Mossa S and Tarjus G 2003 J. Chem. Phys. 119 8069-74 

[6] Karnicky J F, Reamer H H and Pings C J 1974 J. Chem. Phys. 64 4592-600; Kirstein B E 

and Pings C J 1976 J. Chem. Phys. 66 5730-6 
[7] Voronel A, to be published 

[8] Fowler R H and Guggenheim E A 1939 Statistical Thermodynamics (Cambridge) 
[9] Corless R M, Gonnet G H, Hare D E G, Jeffrey D J and Knuth D E 1996 Adv. Comput. Math. 
5 329-59; Hayes B 2005 American Scientist Online 

http : //www. americanscientist . org/template/IssueTOC/ issue/701 

[10] Landau, L D, Lifshitz E M 1980 Statistical Physics (3rd Edition Part 1 ed.) Oxford: Pergamon, 
p. 515 



9 



